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Abstract 

Current fluctuations can provide additional insight into quantum transport in mesoscopic sys- 
tems. The present work is carried out for the fluctuation properties of transport through a pair of 
coupled quantum dots which are connected with ferromagnetic electrodes. Based on an efficient 
particle-number resolved master equation approach, we are concerned with not only fluctuations 
of the total charge and spin currents, but also of each individual spin-dependent component. As a 
result of competition among the spin polarization. Coulomb interaction, and dot-dot tunnel cou- 
pling, rich behaviors are found for the self- and mutual-correlation functions of the spin-dependent 
currents. 
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I. INTRODUCTION 



Noise of transport current through mesoscopic system can provide useful information 



beyond the average current 
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was shown in 



17l |. These studies were focused largely on the total charge current. A relatively 
jject, but now gaining increasing attention, is the spin-dependent current 
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20, |2l| . A number of interesting examples are listed as follows: (i) It 



18l | that the shot noise of spin current is closely related to the spin unit of 
quasiparticles. (ii) In the absence of charge current, the magnetically pumped spin current 
noise through quantum dots in the Coulomb blockade regime can serve as sensitive probe 
to spin decoherence [lOj]. (iii) For a quantum dot strongly coupled to cavity field, shot noise 
of the spin current exhibits clear signature of the discrete nature of the photon states 20|. 
(iv) Spin-dependent shot noise of unpolarized currents can be used to detect attractive or 
repulsive correlations {21]. 

In this work we study the spin-dependent current noises for transport through a pair of 
coupled quantum dots connected with ferromagnetic (FM) electrodes. Owing to the addi- 
tional modulation of the dot-dot coupling, rich noise behaviors are found as a consequence of 
interplay between the Coulomb correlation and the spin polarization of the FM electrodes. 
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lodologically, we employ an efficient particle-number-resolved master equation approach 
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23|, with proper extension by including the spin degrees of freedom. This approach. 



which can be considered as a finite temperature and voltage extension of the "n" -resolved 
quantum Bloch-type equation proposed by Gurvitz and Prager (2^, has the advantages in 
its simplicity of treating properly a broad range of dissipation, as well as its transparency 
in the involving dynamics and many-body interactions. 

The paper is organized as follows. We begin in the next section with the model setup of 
double dots connected with FM electrodes. We then present in Sec. lIIIl a general formalism 
for both the spin- dependent currents and their noises. Numerical results are presented in 
sections HVl and IVl and then followed by the conclusion in Sec.lVIl 
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FIG. 1: Schematic setup for transport through a pair of coupled quantum dots connected with 
ferromagnetic electrodes. 

II. MODEL DESCRIPTION 

The system under study is schematically shown in Fig.[Tl The left and right quantum 
dots are labeled by "L" and "R" , respectively. The total Hamiltonian contains three parts: 

H = Hb + Hs + H'. (1) 

Hb = X]a=L R cr ^afcclfco-^afco- dcscribe the left and right electrodes, which are magnetically 
polarized. Caka (cl^ka-) electron annihilation (creation) operator of the electrode a = 

L or R, with spin a = t or ],. The FM electrodes are characterized by the spin- dependent 
density of states Qatjiuj). Moreover, the densities of the states are assumed to be energy 
independent, i.e. Qaai.^) = Qaa- Then, the FM polarization of the electrodes is simply 
determined by a parameter = {.Qa] — daO/igai + dai), which satisfies —l^Pa^ 1- 
The second Hamiltonian is for the coupled dots: 



a=L,R 

+nJ2idldn. + diM. (2) 

(7 

Here n^o- = dl^^daa and ha = X^o-^qo-) with daa (dl^^) the electron annihilation (creation) 
operator in the quantum dot a and with spin a. Note that here a labels the left and right 
quantum dots rather than the electrodes. Each quantum dot is assumed to have only one 
spin-degenerate energy level -El/r- f^o and U' are the intra-dot and inter-dots Coulomb 
interactions; and Q is the coupling strength of the two dots. The typical value of intra-dot 



interaction Uq is estimated to be 4~5meV, and U' ~ Uo/2 [25 1. 

The third Hamiltonian describes the coupling of the central dots to the electrodes, in 
terms of H' = J2aka {^akCakadaa+ii-C-) ■ Notc that, duc to spiu polarization of the electrodes, 



the tunnel coupling strength would become spin- dependent, i.e. with Ta-\ = (1 + Pa)Ta/'2 
differing from r^j. = (1 — Pa)Ta/2, where = 27r J2k \tak\'^S{£ak — is the total coupling 
strength regardless the spin orientation, and is estimated to be 10^ ~ 10^° s~^ 26|. 

In this work, corresponding to different bias voltages, we will effectively consider the 
following three cases. Case (i): noninteracting (NINT) dots, i.e. each dot can be doubly 
occupied. Case (ii): single-dot Coulomb blockade (SDCB). In this case double occupation 
on the same dot is prohibited, but each of the two dots can hold an electron. Case (iii): 
double-dot Coulomb blockade (DDCB). In this case, there is at most only one electron 
in the two dots. Experimentally, the above mentioned different cases can be achieved by 
appropriately adjusting the applied bias voltage with respective to the intradot and interdot 



charging energies 



23|. 



III. FORMALISM 



In this section we outline the master equation approach for the calculation of the spin- 
dependent shot noise. The main point is to obtain the equation of motion for the conditional 
reduced density matrix of the central quantum dots, which is conditioned not only on the 
transmitted electron numbers, but also on their spins. 



A. Electron-number-resolved master equation 



To achieve the description of spin-dependent transport, we will first extend the "n"- 



resolved quantum master equation [15|, l22|, l23|] to include spin degrees of freedom. Let us 



start with the reduced density matrix, defined as p(t) = TrB[pT(^)], i-e. tracing over the 
electrode states from the total density matrix. Assuming the tunneling Hamiltonian H' is 
weak, the second-order cumulant expansion leads to the master equation 2?! . 

POO 

Pit) = ~iCp{t)~ / dT{c\t)g{t,T)c\r)g\t,T))p{t), (3) 

where the Liouvillian superoperators are defined as CA = [Hs,A], CJ A = [H',A], and 
Q(t,T)A = G(t,T)AG'^ (t,T), with G(t,T) the usual propagator associated with the system 
Hamiltonian Hs- 

To condition the master equation on the electron number transmitted 
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23|, we 



decompose the entire Hilbert space B of the left and right electrodes into the sum of the 



g|^RT' Ri)^ where Uaa is the number of electrons with spin-cr arrived 



subspaces, 

at the electrode a (a=L, or R). Partially tracing over states in each subspace leads to the 
spin-resolved quantum master equation for the reduced density matrix, 



P 



J- "-LT > "RT -J "LT . "RT ■) 

li'^ 4 H n "U ' "Rl -L n "14 ■ "Ri A W rj^ 



,riL^ +1, nr{|, , TiL-j- , np^ , — -l, '^R-f^ 



RL 



_|_ /l( "14. "Rl _|_ /l( )^ 



'"T.I , "HI ~ 1' 



RT ^ ^RL I 



+ H.C. ^ 



Here, -^a(T 



(4) 



E.'Ciliic)d^., and aL^ ^ E.'C'ii(-/:) rf..,, with cSU±/:) = 

J dtC^]^,{t)e^^^^. The involved bath correlation functions are defined as C^]^,(t — r) = 

{fLit)faa'iT)) and C^;i,(t - r) = {faa{t)fl^,{T)), with /c,^ = J2ktc,aCaka and (•■■) = 
TrB[(- ■ ■ )Pb] standing for the usual meaning of thermal bath average. 



B. Spin-dependent currents 



With the knowledge of the above conditional state, the joint probability distribution 

function is obtained as P\("^^'"^^),t\ = Tip where Tr( - ■ ■ ) denotes the trace over the 

system states. It allows us to evaluate the spin-cr dependent current through the junction a 
by 



^ "U ■ "Rl 



LT."RT^ ^ 



eTriV:, 



j-"LT . "RT -J 



where A^^ = Y^^,n^,naj,nj,J^o^^P ""^'""^ ■ Based on Eq. (HD, we obtain 

a' a' 

+[dLp(t)AW-4-Jp(t)4.] + H.c.j. 



(5) 



(6) 



This leads to 



= |eTr{ [dlAl-J - AW4.]pW + H.c.}, 



(7) 
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where p{t) is the unconditional density matrix that satisfies Eq. ([3]), or 

m = -zCp{t) K.' 41p(t) -pW4+]] +H.C.} . (8) 

C. Spin-dependent noises 

Note that the total charge and spin currents through junction a are = Ij^ + I^ and 
J^P = I^ — I^. The total charge current noise and spin current noise can then be expressed 
as 

where the individual spin-dependent noise spectrum is defined as S^'^,{uj) = 
J^jlt cos{ut){{AI^{t),Ai;^,{0)}), with Ai;^{t) = I^{t) - I^. In the following, fiuctuation 



between the same or opposite spin currents is referred to as se^ 

nn 

noise, respectively. Based on the MacDonald's formula [6|, |28|, 



f- or mutual-correlation shot 



29|, it follows that 

5-;(c.) = 2c. dtsmiu;t)-[MZ:it)-&'], (10) 

where is the stationary current obtained from Eq. ([H]), and M^l^t) = 

^ SnLtnLtnR|niy'^ao'^a'o-'-^[(I''^ I'^05 ^] • Usiug the spiu-rcsolvcd mastcr equation (jl]), one finds 

|m::: = le^Tr { [4,4") - A^AA <{t) 

+ [dL4-)+A«4Jp,t5..'5..'+H.c.}, (11) 

where N^{t) can be found from Eq. ([6]), and pst is the stationary solution of the master 
equation (El). Specially, the zero-frequency spin-dependent noise which is of most interest 
reads Q, Q, 19| 



s:i = 2^ [M-:(t) - i:rjt'] i^^. (12) 

So far we have outlined a compact formalism for spin-dependent transport through meso- 
scopic systems. As will show in the following, this approach can be efficiently used to study 
the spin-dependent transport phenomena, say, not only the spin-dependent currents, but 
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TABLE I: Spin-dependent currents, as well as the total charge and spin currents for (i) NINT 
case, (ii) SDCB case, and (iii) DDCB case, respectively. Here, 70 = Fl + Fr, 7s = 2Ti^ + Fr, and 
7(j = 4Fl + Fr are the corresponding total effective tunneling rates for these three cases. 

Cases (i) (ii) (iii) 

n/p n I r,\ rLrRnV7o i+p 2rLrRnV7, i+p 4rLrR.nV7d 
/e l-L^/^;f^2+i(i+p)2rLrH, 2 n2+i(i+p2rR,/7s)rLrH, 2 2c5T(TTF){fZr^77^ 

ri/p (^ ^i^^n^Vyo i-p 2rLrRnV7s i-p 4rLrRnV7d 

J /e li /'Jf^2+i(i_p)2rLrR 2 n2+i(i+p2rj,/7,)rLrR 2 2fPT(T+OTZr577^ 

fch/p 8rLrRo2[(i-p^)rLrR+4o2]/7o 2rLrRnV7, 4rLrRnV7d 

' [{i+p2)rLrR+4n2]2-4p2r2r2^ n2+i(i+p2rj,/^,)rLrR 2n5T(T+OTZr|7My 

fsp/„ 8rLrRn2[(p2„i)rLrR+4n2]/7o 2rLrRn2/7, 4rLrR!^2/^^ 

/ ^ [(i+p2)rLrR+4Q2]2_4p2r2r2 P^2_^i(^_^p2p^^^^^j,^p^ ^^^WT(T+^iW^jfhI) 

also their noise properties. The major approximation involved in our formalism is the stan- 
dard second-order Born approximation with respect to the perturbative expansion of the 
tunneling Hamiltonian. Apparently, this approach is applicable under the Markov-type con- 
dition for relatively high temperatures, say, higher than the tunneling width [s^. In this 
case, the formalism can be safely applied for arbitrary bias voltages, including those with 
Fermi level of electrode in near resonance with the energy level of the central (transport) 
system. Also, the present formalism is applicable to low temperature under large bias volt- 
age, i.e. the electrode's Fermi surface is further away from the central system's energy level 
than the level broadening. Similar statement holds also for the widely applied rate equation 
approach at zero temperature [24 1. Finally, at low temperature and under low bias voltage 



(near-resonance transport), an improved self-consistent Born approximation 3l| or exact 
hierarchical equations of motion approach [32i] can be satisfactorily employed. 

IV. CURRENT CHARACTERISTICS 

For simplicity, we assume the left and right electrodes have the same spin polarization 
Pl = Pr = P, and the two levels are in resonance, i.e. Ei^ = E^. This can be achieved, 
in experiment, by applying appropriately gate voltages on the left and right dots (see, for 



instance, Ref. 33|). Hereafter, we consider large bias voltage and low temperature in each 
case, which makes the Fermi functions relevant to the transport processes be either one or 
zero. Under this condition, electrons only transport in one direction, say, from left electrode 
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Polarization p Polarization p Polarization p Polarization p 



FIG. 2: Currents versus spin polarization p for (i) NINT case (solid curves), (ii) SDCB case (dashed 
curves), and (iii) DDCB case (dotted curves), respectively. Figures (a)-(d) show respectively, the 
spin-up, spin-down, total charge and total spin currents for $7 = 0.2r, and those in (e)-(h) are the 
corresponding currents for 17 = 2. OF. 

to the right one. 

Listed in table [T] are the stationary spin-dependent currents, together with the total 
charge and spin currents. Numerical results are presented in Fig. [2], where we further as- 
sume symmetric coupling, i.e. Fl = Fr = F. Here, the currents are plotted against the spin 
polarization under the condition of large bias voltage. At finite bias, the energy renormaliza- 
;ion will be involved in the present formalism. It arises from the imaginary part of the rate 
3o| , and can be included readily within our theory by following 34 , 35 1 . The main features 
found here are summarized as follows: (i) the currents are suppressed in the presence of 
Coulomb interactions; (ii) there exists a symmetry relation between V and /-'-, which makes 
us only need to consider one of them, e.g. V] (iii) the spin-up currents {V) and the total 
spin currents increase monotonically with p, while the charge currents reach their maxima 
at p = 0. 

The results for the noninteracting and weakly coupled dots are shown by the solid curves 
in Fig.[2](a)-(d), which differ from the behaviors stated above in (iii). In this weak dot-dot 
coupling regime, the resultant large occupation probability of the left dot leads to simple 
expression for the currents. For instance, the spin-up current in this case can be well 
approximated by P ~ TLfPLi with plo- the probability of the left dot being occupied by a 
spin-cr electron. As the polarization of the FM electrodes is gradually altered from p = —1 
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FIG. 3: Individual spin-dependent noise for NINT case (solid curves), SDCB case (dashed curves), 
and DDCB case (dotted curves), respectively. Figures (a)-(c) are the noises for weak coupling 
(0 = 0.2r), and those for strong coupling [0. = 2. OF) are shown in (d)-(f). 

io p = 1, decreases rapidly, while Fl^ increases linearly. Then, a turnover behavior for 
the spin-up current is formed as shown by the solid curve in Fig.[2](a). This is also the basic 
reason for the unique features of the charge and spin currents shown by the solid curves in 
Fig.EKc) and (d). 



V. NOISE CHARACTERISTICS 



23| . we will investigate various 



In this section, instead of using the circuit noise 
fluctuations of spin currents through left or right junction, i.e. the auto-correlations 5"; 
with a = L or R. For the cross-correlations, they simply satisfy the relation = —S] 



(T(T 

aa 



in the present two-terminal case, as we have checked. It should be noted that for a three- 
terminal structure such a simple relation generally does not take place (see, for instances. 



Ref. 



36 



37 



38|). 



A. Self-correlation shot noise 



We now turn to the noise characteristics, and first to the self-correlation shot noise. Note 
that there is a symmetry between S]J^ and S"^^ with respect to the spin polarization, i.e.. 
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Fig.[3](a) versus (b) and (d) versus (e). Therefore, it needs only to consider either S^J^ or 
S^l^. For SjJ^ in Fig.[3]^a) and (d), it is found that both the Coulomb correlation in the dots 
and the spin polarization of the electrodes will enhance the Fano factor. The shot noise is 
of NINT < SDCB < DDCB, consistent with the Coulomb correlation strength. The only 
exception in terms of spin polarization is that shown by the solid curve in Fig.[3]^a), for the 
weak dot-dot tunnel coupling and in the absence of many-body Coulomb interaction. In 
this case the Fano factor does not monotonically increase with the polarization degree of the 
electrodes. This behavior originates from the same reason as that leading to the shoulder 
behavior of the charge current in Fig.[2](c). 

Remarkably, a profound, strong super-Poisson behavior can be developed by increasing 
the dot-dot tunnel coupling for both the SDCB and DDCB cases, provided the electrodes are 
properly spin polarized, see the dashed and dotted curves in Fig.[3](d). This novel behavior 
can be understood in terms of the so-called dynamical spin blockade mechanism. Take Sj^^ 
again for illustration. If the electrode is sufficiently spin-up polarized, the spin-up electrons 
can more easily pass through the two dots than the spin-down electrons. Furthermore, 
the SDCB (DDCB) does not allow double occupation in single dot (double dots). Then a 
mechanism of fast-to-slow channels is developed, which results in a bunching behavior and 
also the super-Poisson statistics of the transport electrons. 

Note that the strong Coulomb interaction is essential to the bunching behavior of tun- 
neling events, otherwise the spin-up and spin-down electrons will transport independently, 
and will not result in the super-Poisson statistics at all [1]. Also, a relatively strong dot-dot 
coupling is required for the super-poisson noise in the coupled double dots system. If the 
dot-dot coupling is weak, electron will attempt to stay longer in the left dot. This will 
weaken the bunching behavior, thus lead all the noise components to sub-Poisson, see the 
dashed and dotted curves in Fig.[3](a) and (b). 

The mechanism of dynamic spin blockade plays also an essential role in a three-terminal 
quantum dot [371, l38| , and in an interacting quantum dot with intradot spin-flip scattering 



connected to FM electrodes 



391]. Meanwhile, a similar mechanism called dynamic charge 



blockade is responsible for the super-Poisson noise in interacting two-channel systems 



41 



45|, as well as double quantum dots structures |42l. |43|. 
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B. Mutual spin correlation shot noise 



We now turn to the mutual spin correlation shot noise {S]J^ or S]J^), which are symmetric 
to the polarization, as shown in Fig. [31(c) and (f). For noninteracting dots, the mutual 
correlation between opposite spin currents are zero, as shown by the solid curves in Fig.[3](c) 
and (f). This is simply because the spin-up and spin-down currents are uncorrelated in the 
absence of Coulomb interaction, despite the electrodes are spin polarized. Similar behavior 
was also found in where the mutual correlation is zero for noninteracting single dot 
connected with normal electrodes. 

For interacting dots, i.e., the SDCB and DDCB shown in Fig.O (c) and (f), the mutual 
correlation of spin-up and spin-down electrons reveals relatively complex features which can 
be either positive or negative. The main features are as follows, (i) For SDCB, negative 
mutual-correlation is found for the weakly coupled dots, if the electrodes are weakly or 
moderately polarized; but positive correlation can be formed if the electrodes are strongly 
polarized, (ii) Also for SDCB, the mutual correlation is fully negative for strongly coupled 
dots, independent of the polarization degree of the electrodes, (iii) For DDCB, in the case 
of weak dot-dot coupling, the mutual correlation is fully negative, (iv) Again for DDCB, 
but in the case of strong dot-dot coupling, the mutual correlation is positive and can be of 
strong super-Poisson for electrodes sufficiently polarized, while it is negative for electrodes 
weakly or non-polarized. 

In this context, we notice that the super-Poissonian self-correlation does not necessarily 
imply a positive mutual correlation, as shown by the dashed curve in Fig.[3]^f ) for the SDCB, 
where the mutual correlation, despite the polarization degree of the electrodes, is fully 
negative. A simple and unambiguous understanding for the above positive and negative 
mutual correlation is yet to be found. It seems more subtle than the super-Poisson statistics 
of the self-correlation of current. 

C. Total charge and spin current noises 

The total charge and spin current noises are nothing but the combination of those compo- 
nents of the spin-dependent noises, according to Eq. The results are displayed in Fig.Hl 
where most features such as the sub-Poisson and super-Poisson behaviors can be accordingly 
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TABLE II: Spin-dependent noise, as well as the total charge and spin current noises for p = in 
differ ent cases. 

Noises S]Ja {Saa) Saa {Saa) ^aa Saci 

case (ij eJo— (f5+4jpp— U ^^^o — (f^+if^^p— ^eio — (f3+4fpp— 

pflsp (\\^ pT 9r^+i4r"n2+28n4 (iir2+4n2)n2 9r4-8r^n"+2oa^ 

rasp niil pLrl±6r!n!+84nl 2pL(IE^±5^^ grL r^-sr^n^+esn^ 

case (mj eid (r^+ion^)^ ~^^^d-(f3+xo!Fp~ -^eid (r^i+ion^)^ ^^■'d 



understood in terms of the above interpretation to the partial noises. 

A noticeable result is the spin current noise resulting from the DDCB. We find that it is 
constantly Poissonian, regardless the polarization degree and dot-dot coupling, as shown by 
the dotted curves in Fig.Hl^b) and (d). Similar Poissonian spin current noise was also found 



m 



2l| . for single quantum dot in the Coulomb blockade regime. 



To complete this section, we consider all kinds of noises discussed above in the limit of 
p = 0. In this unpolarized situation, simple analytic results can be obtained, as shown in 



12|, 



[Sland Fig.H The total 
4J| by an overall factor 



table Hn which serves as a complement to the numerical results in 
charge current noise in the NINT case differs from the ones in [l^. 
of 2 (implied in /q, see table H]). This factor originates from the spin degree of freedom, which 
was ignored there. For interacting cases, under symmetric coupling (Fl = Fr) and resonant 
levels (i?L = -Er), the total charge current noises are found here to be sub- Poissonian, 
in either the SDCB or the DDCB cases. Very interestingly, however, if Ei^ ^ E^i and 
Tl > Tr, the interplay of quantum coherence and the DDCB would develop a dynamical 
charge blockade mechanism, which then result in a remarkable super-Poissonian noise. Our 
result in the DDCB case slightly differs from the ones in 10|, |42]- Again, the difference arises 



from the spin degree of freedom, which was neglected in the previous work. 



VI. CONCLUSION 



To summarize, based on an efficient particle-number resolved quantum master equation 
with inclusion of the spin degrees of freedom, we investigated the spin-dependent noises in 
transport through a pair of coupled quantum dots connected with ferromagnetic electrodes. 
The modulation of the dot-dot coupling, and the interplay between Coulomb interactions and 
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FIG. 4: Total charge and spin current noises, obtained by appropriate combination of the individual 
noise components as shown in Fig. [3] according to Eq. ([9]). The plot parameters are the same as in 
Fig.El 

spin polarization in the electrodes give rise to rich noise behaviors, such as super- Poissonian, 
constant Poissonian noises, as well as positive and negative mutual correlations. These 
unique noise features can serve as additional tool in experiments for revealing the intrinsic 
dot-dot coupling, as well as the Coulomb interactions that involved. The results presented in 
this work were carried out at zero temperature. At finite temperatures, backward processes 
are possible, which will in general reduce the current, but enhance the noise, as clearly 
demonstrated in 



45|. 



Finally, we briefly discuss the measurement of spin-dependent noise. It was shown in [46 1 
that in hybrid ferromagnetic-normal metal structures, the spin current noise exerts a fluc- 
tuating spin torque on the magnetization vector, which causes an observable magnetization 
noise. By measuring this magnetization noise, which has been realized experimentally 
one then obtains the spin current noise. 
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